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We point out that the lower bound on the primordial magnetic field required to seed the galactic 
dynamo is significantly relaxed in an open universe or in a universe with a positive cosmological 
constant. It is shown that, for reasonable cosmological parameters, primordial seed fields of strength 
10""^" Gauss or less at the time of galaxy formation could explain observed galactic magnetic fields. 
As a consequence, mechanisms of primordial magnetic seed-field generation that have previously 
been ruled out could well be viable. We also comment on the implications of the observation of 
micro- Gauss magnetic fields in galaxies at high redshift. 
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Magnetic fields pervade most astrophysical systems [0 , 
but their origin is unknown. Spiral galaxies are observed 
to possess large-scale magnetic fields with strength of the 
order of 10~^ G and direction aligned with the rotational 
motion. A plausible explanation is that galactic mag- 
netic fields result from the exponential amplification of 
an initially weak seed field by a mean- field dynamo . 
Many proposals have been put forward regarding the ori- 
gin of such a seed field. One suggestion is that it might 
arise spontaneously from non-parallel gradients of pres- 
sure and charge-density during galaxy formation Q. A 
wider range of possibilities is offered if the seed field is 
of primordial origin. This category includes cosmological 
magnetic fields Q as well as magnetic fields created by 
any of a number of early-universe particle-physics mech- 
anisms such as collisions of bubbles in a first-order 
phase transition Q or false- vacuum inflation 

The seed-field strength required at the time of com- 
pleted galaxy formation (igt) for a galactic dynamo to 
produce the present magnetic field strength ~ 10~^ G 
is usually quoted in the range ~ 10~^^ - 10~^^ G. Such 
lower bounds are obtained by considering the dynamo 
amplification in a flat universe with zero cosmological 
constant for "typical" values of the parameters of the 
aw-dynamo. The seed field must also be coherent on a 
scale at least as large as the size of the largest turbulent 
eddy, ~ 100 pc [1). Most proposed models of primordial 
seed-field generation fail to meet these requirements as 
formulated above. 

In this paper we address the issue in the light of re- 
cent developments in cosmology. Observations of distant 
type-IA supernovae and of anisotropics in the cosmic 
microwave background (CMB) ||ll[] in combination have 
made it increasingly likely that the universe is less dense 
than the critical density and has a positive cosmological 
constant A. Most previous studies of magnetic fields have 
assumed a A = universe with critical matter density. 
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We shall recalculate the constraints on primordial seed 
fields for general Friedmann universes with matter den- 
sity parameter JIq and vacuum energy density parameter 
Ao = A/{3Hq) such that f^o + -^o < 1 (the subscript 
here indicates quantities at present time and Hq is the 
Hubble parameter) . In addition to finding revised bounds 
on the seed field at time tgf we shall trace the evolution 
of the magnetic field back to the time of radiation de- 
coupling, tdcc- Prior to decoupling the evolution of the 
magnetic field proceeds via complicated plasma processes 
and depends on the field's initial strength and correlation 
length After decoupling there is sufficient residual 

ionisation for the magnetic field to be frozen into the 
plasma; the evolution is simple and independent of the 
mechanism of generation. Thus idee is a natural epoch 
for imposing bounds on primordial magnetic fields. 

We begin by considering the aw-dynamo ||^,^ which is 
a well-studied model of amplification of magnetic fields. 
It is powered by the differential rotation of the galaxy 
in combination with the small-scale turbulent motion of 
ionised gas. By separating the magnetic field into a large- 
scale mean field B and a random, turbulent field b one 
obtains a system of equations with exponentially growing 



solutions Bi, 



for the azimuthal component of 



S Electronic address: A.C.Davis@daintp.cain.ac.uk 
Electronic address: M.J.Lilley@daintp.cain.ac.uk 
Electronic address: 0.Tornkvist@damtp.cam.ac.uk 



the mean field in the plane of the disc. The dynamo 
amplification rate F appears as an eigenvalue that must 
be determined numerically. 

Unfortunately the value of F is rather sensitive to the 
parameters of the dynamo model |]l3| , which include the 
root-mean-square velocity and magnetic diffusivity of the 
turbulent plasma as well as the angular- velocity gradient 
r duj{r)/dr. For reasonable estimates of these quantities, 
one obtains F^^ in the range 0.2 < F~^ < 0.8 [Gyrs]. 
Because of the exponential growth, such an uncertainty 
quickly translates into an uncertainty of many orders of 
magnitude in the total amplification, which may or may 
not rule out various seed- field mechanisms. The point 
of this paper is not to linger on these uncertainties but 
rather to emphasise the tremendous increase in amplifi- 
cation that will occur in an open universe or in one with 
a positive cosmological constant for any value of F. We 
shall present results for the two values that appear most 
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frequently in the literature, F -'^ = 0.3 Gyrs and 
r-i = 0.5 Gyrs §. 

Any magnetic seed field is exponentially amplified until 
it reaches the equipartition energy (B ^ few /iG) when 
further growth is suppressed by dynamical back reaction 
of the Maxwell stresses on the turbulence. Assuming 
that the dynamo mechanism begins to operate around 
the time of completed galaxy formation igf, the lower 
bound on the strength of the seed field at this epoch is 
given by 



Bgf > Boe 



-r(to-tgf) 



(1) 



where Iq is the age of the universe. Expressions for to in 
different cosmologies are found in e.g. Refs. and [p^ . 
In particular, for a given value of Hq, open universes and 
universes with a positive cosmological constant are signif- 
icantly older than the Q = 1 Einstein-de Sitter universe. 

The time of galaxy formation t^t can be estimated from 
a spherical collapse model. As we shall show presently, 
galaxies of a given average density Pgai have collapsed 
at approximately the same time after the Big Bang for 
all realistic cosmological models. Galaxies in an open or 
A > universe are therefore older, giving the dynamo 
mechanism more time to — tg[ to operate. Consequently, 
a much smaller magnetic field Bg{ can seed the dynamo 
and still give the observed micro-Gauss field Bq. 

The spherical collapse model describes the non- 

linear collapse of a bounded spherical region with average 
local density pi larger than the critical density at some 
initial time ti in the matter-dominated era. This over- 
density causes the sphere to break away from the Hubble 
expansion, reach a maximal (turn- around) radius r^, and 
eventually collapse to form a gravitationally-bound sys- 
tem. The general equation of motion for a shell of radius 
r enclosing mass M is 



GM 



-Ar^ = E , 
6 



(2) 



where E is a constant. This equation can be separated 
to yield 



3/2 



dt = ± 



V2GM ^1 - (1 + f3)x + I3x^ 



(3) 



where x = r/r,„ and j3 = Arj^/(6GAf). The exact so- 
lution can be expressed in terms of incomplete elliptic 
integrals , but we choose instead to expand the right- 
hand side of Eq. in (3, obtaining the more convenient 
parametric solution 



where T is a constant which can be neg lected [0- We 
see that turn-around occurs at a time tm correspond- 
ing io 6 — IT. As the spherical region recollapses for 
t > tiji, random non-radial particle velocities become 
important; the simple collapse model breaks down and 
the collapse is halted at a final radius rvir given by the 
virial theorem. For universes with zero cosmological con- 
stant Tvir = "Tin/S. For < /3 < 1/2 (which is re- 
quired for collapse to occur) Lahav et al. 15 1 showed 
that 1/2 > Tvir/?"!!! > 0.366 and also obtained the ap- 
proximate relation 



^vir 



1- /3 

2- /3 



(6) 



We can estimate (3 for a galaxy: Taking < Aq < 1, 
M = 1O"M0 w 2 X lO^^g and rvir ^ 15 kpc we get 



2GM (0.366)^ 



' 3 X 10" 



(7) 



The small value of f3 signifies that the vacuum energy 
density plays a negligible role compared to the matter 
density in the collapse of objects as small and dense 
as galaxies. From Eq. (^ it follows that we can set 
Tvii. = for all realistic values of the cosmological con- 
stant. Moreover, we can neglect the /3-dependent terms 
in Eq. (1). 

It is generally assumed |l^ that gravitational collapse 
is complete at the time tvh- > im when r approaches zero 
in Eq. (^), corresponding to sa 27r. |^ This assumption 
is supported by N-body simulations and, because of the 
small value of (3, remains valid in any realistic Friedmann 
cosmology. From Eq. (||) we then have ^vir ~ "^tm as well 
as tl^ = 37r/(32Gpm), where = 3M/(47rrj^) is the 
average density of the spherical region at turn-around. 
It follows from rvir — ''m/2 that pgai = Spm. Finally, 
with tg{ — tvir, all these relations combine to give 



Pgal 



37r 



(8) 



This is the result that we have sought. It shows that 
the relationship between the average density of a galaxy 
Pgai and the time of completed galaxy formation tg{ is in- 
dependent of cosmology. The current galactic density is a 
quantity which can be measured using methods that also 
do not depend on cosmology. The redshift of galaxy for- 
mation Zg{ corresponding to a given value of pg^i is given 
implicitly by well-known time-redshift relations (see e.g. 



t + T = 



y-(l-COS( 

3/2 



VGM 



- sin e) 



— (66 e - 93 sin 61 -f 15 sin 261 - sin 361) + 0(f3^) 
96 



(4) 



(5) 



^The naive estimate, that collapse is complete when the ra- 
dius r in Eq. (^ reaches the virial radius rvir (corresponding 
io9 = 3n/2), is unrealistic as the radius decreases more slowly 
during virialisation than the spherical collapse model would 
imply. 
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FIG. 1. Lower bound on the seed field at galaxy formation Bgf vs Q,o: (a) universe with A = 0, (b) flat A universe. 



Refs. |TJ] and |jT^). It may at first seem mysterious that 
ilo does not enter in Eq. ^ or any of the derivations 
leading to it. The reason is that (for A = 0) the same 
average local density (larger than the critical density) is 
required for a spherical region to collapse regardless of 
the density of the surrounding universe. By Birkhoff's 
theorem the evolution therefore proceeds in an identical 
manner. 

We are now in a position to calculate bounds on mag- 
netic seed fields Bgi at the time of completed galaxy for- 
mation in different cosmologies. We take Bq = 10^^ G 
and pgai = 10^^'' gcm"^. The latter value corresponds 
to the average density of the galactic halo rather than 
the central disc, whose density is ~ 10~^^ gcm~^. The 
reason for this choice is that the spherical collapse model 
uses the simplified assumption of a uniform-density "top- 
hat" profile of the galactic density distribution, and since 
the halo comprises most of the volume of the galaxy, this 
value seems more appropriate. The precise value of pgai 
is of little importance as our results are quite insensitive 
to it. 

The results are displayed in Fig. |](a) for a A = uni- 
verse and in Fig. ^(b) for a flat A universe {Qq + Xq = 1). 
The quantity h is the Hubble parameter Hq in units of 
100 km sec~^ Mpc""'^. For comparison, the straight hor- 
izontal line in each plot shows the constraint of Bg{ > 
10~^° G given by Rees It can be seen on these 

graphs that in an open universe, and particularly in a 
universe with a significant cosmological constant A, this 
requirement is too strong. For reasonable cosmological 
parameters and the same value of F, the dynamo mecha- 
nism could generate currently observed galactic magnetic 
fields from a seed field of the order of 10"'^" G or less 
at the completion of galaxy formation provided that the 
seed field is coherent on a scale ^gf ^ 100 pc. 

We shall now evolve these bounds back to the time 
of radiation decoupling, taking the conservative view 
that there is no magnetohydrodynamic turbulence or 



dynamo mechanism operating during gravitational col- 
lapse (see, however, Ref. for more optimistic propos- 
als). The magnetic field is assumed to be frozen into the 
plasma and its evolution is determined by flux conser- 
vation Br^ = const., where r is a length scale evolving 
with the matter, i.e. r ~ (p)~^/^. Care must be taken 
not to associate this length scale with the scale factor 
a{t), as a collapsing galaxy is decoupled from the Hubble 
expansion. One obtains B 



B, 



gf 



B, 



doc 



Pgal 
Pdcc 



2/3 



2/3 



1 



Pgal \ 

PO J (1 + ZdGc)^ 



(9) 



where we have used the energy conservation relation 
p — /3o(l + 2)"^ for the matter component and the fact 
that the matter density pdoc at the epoch tdoc is very 
nearly uniform. The redshift of radiation decoupling, 
Zdec, has a weak dependence on the product in- 
volving the fraction of critical density in baryons b and 
is constrained to lie in the interval 1200 Zdoc ^5 1400 

Note that the magnetic field will decrease between tdoc 
and <gf , since by virtue of the Hubble expansion the phys- 
ical volume of the galaxy is larger than the volume con- 
taining the same mass at ^doc- The depletion depends 
on the cosmological parameters via the present matter 
density po ~ 1-88 x lO'^^f^o^^ [g cm'^]. It can be 
seen that the depletion is somewhat smaller in universes 
with ilo < 1- This further increases the effect of cosmo- 
logical parameters in relaxing the bounds on primordial 
seed fields. The resulting bounds for B^cc are shown in 
Fig. ||(a) and Fig. ||(b) for a A = universe and for a flat 
A universe, respectively. 

We shall now address the issue of the correlation length 
of the magnetic field. In order for the galactic dynamo to 
begin to operate, the correlation length of the seed field 
at the time of completed galaxy formation must satisfy 
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Cgf ~ 100 pc 1^. 1^ Using the spherical collapse model, 
one can calculate the physical scale r^oc at the time of 
radiation decoupling that will evolve into the size of a 
galaxy. At any time before the onset of gravitational col- 
lapse the matter density follows the Hubble expansion 
and it makes sense to express rdcc in the constant co- 
moving quantity x defined by r = a{t)x. The comoving 
scale X corresponding to a galaxy is given by 



Xgai = 0.95 {VLoh^) 



[Mpc] 



(10) 



where M12 = M/IO^'^Mq. 

The correlation length ^ can be written as a fraction 
of the radius of the galaxy, ^ = jyrvir. With the sim- 
plified assumption of the spherical collapse model that 
the collapsing region has uniform density, the collapse 
is homogeneous and isotropic and different scales col- 
lapse proportionately. Assuming that the magnetic field 
is frozen into the plasma between tjec and t^f we have 
a^corr = f?a;gai. For a galaxy, M12 « 0.1, and the typi- 
cal length scale of the turbulent motion, ftmb = 100 pc, 
corresponds to 77 « 1/150, giving the following bound on 
the comoving correlation length 



a;turb = 5-10 kpc 



(11) 



for observationally realistic values 0.25 > Q,Qh^ > 0.025. 
This bound is somewhat higher than that stated in 
Ref. iQ. The bound should not be applied before idee, 
since the correlation length then evolves according to 
complicated magnetohydrodynamic processes and is not 
proportional to the scale factor a{t) [|l2|. 



A more conservative bound, used by many authors, is 
^gf > 1 kpc. 



In general, primordial seed fields produced by particle- 
physics or field-theory mechanisms are too incoherent to 
meet the requirement posed by Eq. (1^). However, there 
is a possibility even for a less correlated magnetic field 
to pass the requirement provided that it has sufficient 
strength to satisfy the bound on B^cc after r.m.s. coarse- 
graining over the scale given by xtmh ■ The said procedure 
results in an root-mean-square field 



Bit 



■^dcc 



d/2 



B. 



dec ; 



(12) 



where iJrms is the quantity that must satisfy the bound 
given in Fig. with i?dec and Xdec being the strength 
and comoving correlation length, respectively, of the pri- 
mordial seed field evolved from formation to idee- The 
exponent d can equal 1, 2, or 3 depending on the aver- 
aging procedure used. This complicated issue (20| shall 
not be addressed in this paper. 

There have been observations of micro-Gauss fields at 
redshifts of z = 0.395 and z = 2 ||2|, ahhough the 
latter has been criticised ||2^. If correct, these observa- 
tions are difficult to explain in a flat universe with A = 0. 
They may, however, be easier to understand in an open 
or A universe. Applying our model to the z = 0.395 case, 
with -Bo. 395 = 10~^ G, we obtain for a flat A universe the 
bounds at t = idee shown in Fig. ^. Hence a seed field 
of 10~^° G at tdoc, or equivalently lO^-^'^ G at igf, could 
account for this observation. 

If we attempt a similar analysis in the z = 2 case, the 
required seed field is sufficiently high that it would have 
other cosmological irnplications, e.g. on the CMB j24j or 
structure formation Consequently we conclude that, 
unless the dynamo parameters are radically different for 
high column density Ly-a clouds (e.g. if they have fast- 
spinning cores and thereby have a higher angular- velocity 
gradient \rduj/dr\ [p6l) these observations cannot be ex- 
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plained by amplification of a primordial seed field by a 
galactic dynamo. 

In this paper we have reconsidered the constraints on 
the primordial magnetic field required to seed the galac- 
tic dynamo in the light of recent cosmological advances. 
We have shown that, in an open universe or a universe 
with A > 0, a much smaller seed field is required to ex- 
plain the observed micro-Gauss fields in galaxies. As a 
consequence, mechanisms of primordial magnetic seed- 
field generation that had previously been ruled out, on 
the grounds of giving too small strength or correlation 
length, could well be viable. We have evolved the bounds 
back to the epoch of radiation decoupling idee, assuming 
that from tdcc to the present the magnetic field is frozen 
into the plasma and evolves first via flux conservation 
and thereafter by amplification via a dynamo mechanism. 
The remaining problem is to evolve primordial magnetic 
fields from the time of their generation to idee taking 
into account various plasma effects. A step in this direc- 
tion has been taken in Ref. This work needs to be 
generalised to different cosmologies, although it can be 
expected that the main cosmological effects occur at late 
times. 
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